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Abstract 

The exact expression of the free energy of a planar array of a Ginzburg- 
Landau chains with nn and nnn interaction is obtained. The critical 
behaviour of the specific heat is not qualitatively modified by taking into 
account the nnn interaction. 

The Ginzburg-Landau functional, sometimes called Ginzburg-Landau- Wilson 
effective Hamiltonian, is of utmost interest in physics. In 2D, the real Ginzburg- 
Landau fields may describe uniaxial ferromagnets [T] or ferroelectrics [3] , ultra- 
thin films of hydrogen-bonded chains [3], [4] etc. Recently, the free energy of a 
planar array of Ginzburg-Landau chains, with near-neighbor (nn) interaction, 
has been calculated exactly and a 2D Ising-like critical behaviour of the specific 
heat has been obtained [5]. This result has been extended to the 3D case in [6] 

In this paper we shall extend the approach of [5], taking into account the 
next-near neighbor (nnn) interaction between the chains of a planar array. The 
Ginzburg-Landau functional of the system is: 
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where the field ipj(x), on the j— th chain, is real, and satisfies cyclic bound- 
ary conditions, ipj(x) = ipj+]\f(x). The parameters a, 6, c, £o have their usual 
signification, see [5] for details. 

In order to obtain the thermodynamics of our system ([T]), we shall follow, 
as we did in [5], the approach of Scalapino, Sears and Ferell [7]. According to 
these authors, the free energy of H} is proportional to the ground state energy 
of a Hamiltonian associated to ([T]) - the so-called transfer matrix Hamiltonian. 
In our case, the transfer matrix Hamiltonian, associated to |T]), is: 
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It describes a ID system of anharmonic quantum oscillators, coupled by quadratic 
nn and nnn interactions; so, it presents also an intrinsic physical interest [3]. 
The "mass" in ^ is temperature-dependent: 
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It is convenient to replace ^ with a simpler Hamiltonian: 

Htm — k B TT n HTM 



where: 
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Tmf is the "mean-field transition temperature", which, in fact, does not corre- 
spond to any physical phase transition, but is just one of the parameters of the 
Ginzburg-Landau theory, describing the temperature dependence of a : 



In fact, 
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The Fourier transform of the interaction potential, 
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where clq is the lattice constant of the chain, plays an important role in this 
approach, because the ground state of (0 can be written as (see [5] for details): 
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where J7o is the ground state of the one-particle anharmonic oscillator, 

£7o is T— dependent, as a function of A, according to eq. JHJ). Finally, the 
ground state energy of ([7]) can be written as: 



£ G = gfe ,/l + 2 ^ + C 2 2) sin 2 tf) (1 - fcf sin 2 (13) 
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where fc 2 , fc| are parameters depending on c\ , c 2 and f2o . If 

ci > 0,c 2 < 0,ci > |c 2 | (14) 

they satisfy the conditions: 

fc 2 > k\ > 0; fc 2 -> i/ c 2 -> (15) 
So, the free energy of our system (pj is proportional to the integral: 

E[k u k 2 ) = Jj ^{1 - k 2 sin 2 9) (l - k 2 sin 2 6) d9 (16) 
More precisely, 



F = Nn k ^Jl + 2 -^±^E(k 1 ,k 2 ) (17) 
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The function E(ki,k 2 ) can be considered a generalization of the complete 
integral of second kind, E(k). For fc 2 = 0, corresponding to the case when the 
nnn interaction is neglected, E (fcj.,0) = E(k\), and (fTT|) reduces to our previous 
result, eq. (69) of [5]. 

Using (252.19) and (262.17) of 0, E(ki,k 2 ) can be put in the form: 

£ ( fc "*'>°- 2 tt 4-l)VlW <18) 
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where 



a* = < o fc 2 = ^ ^ (19) 
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and the conventions of [5] for elliptic integrals and functions have been adopted. 
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The formulae (fTT)) . (fTS|) give the exact expression of the free energy of our 
system, ([1]). The addition of an external "magnetic" field H can be made very 
simply, adding a linear term in the anharmonic Hamiltonian (|12[k in this way, 
flo will become field-dependent, Qq(T,H). This is an important advantage of 
this simpler model, compared to 2D Ising, where the addition of an external 
field defies an exact solution. 

Let us shortly discuss now the critical behaviour of the system. According 
to (fTT|) , (fT5|) , the free energy is proportional to: 



fc 2 (1 - fef) E(k) + k\ (l - kl) K{k) + {k\k\ -k\- kf) n(a 2 , k) (20) 



The singularities of this expression might appear at k = 1, due to the terms 
in K{k) and II(a 2 , k). In fact, the only singularity of II(a 2 , k) (see [8] (904.00)) 
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is due to the K term, because Heuman's lambda function can be written as ([8 
(904.01)) 
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where 



a =E, a 2 \ (2K-E)k 12 , a 4 = I (AK - 3E) fc' 4 , ... (23) 
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and tends to 

-(^ if fc 2 -> 1 or fc' 2 = 1 - fc 2 -> (25) 

So, it is easy to see that, in (|2H)) . for fc 2 — > 1, the AT terms compensate each 
other, and the most singular contribution, coming from E and from the second 
term in (|16p . is proportional to 

fc' 2 lnl (26) 

The cancellation of the most singular terms in (|2"D|) is quite similar to that 
obtained by Fan and Wu in their calculation of the specific heat of a 2D Ising 
model with nnn interaction. The term (|26|) gives a logarithmic singularity of 
the specific heat, which occurs at a critical temperature given by the equation 
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k = 1 or, equivalently (see eq. (|19j) ). k\ = 1. So, the presence of a nnn term 
does not change qualitatively the critical behaviour of the system. 

The influence of the nnn interaction on the critical behaviour was subject 
of intense debate for the Ising model (as it is quite generally accepted, the 
planar array of Ginzburg-Landau chains belongs to the Ising universality class). 
The conclusion of several analytical approximations ([S], [TUj) and Monte Carlo 
simulations ([11], [12]) is that the nn Ising critical behaviour is not qualitatively 
modified by the nnn interaction, at least for small values of the coupling constant 

C2- 

Our result agrees with these conclusions. As §5$ is also the ground state en- 
ergy of a chain of quantum anharmonic oscillators with nn and nnn interaction, 
it can be also used in the study of quantum phase transitions in this system, 
which, apart of being a transfer matrix Hamiltonian, presents a certain physical 
interest in itself. 
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